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Abstract 

r*j | We study a system consisting of a heavy quantum particle, called tracer particle, coupled 

to an ideal gas of light Bose particles, the ratio of masses of the tracer particle and a gas par- 
ticle being proportional to the gas density All particles have non-relativistic kinematics. The 
tracer particle is driven by an external potential and couples to the gas particles through a pair 
potential. We compare the quantum dynamics of this system to an effective dynamics given 
by a Newtonian equation of motion for the tracer particle coupled to a classical wave equation 
^ . for the Bose gas. We quantify the closeness of these two dynamics as the mean-field limit is 

£> \ approached (gas density — » oo). Our estimates allow us to interchange the thermodynamic with 

the mean-field limit. 

o; 

o 1 Introduction 



As a model for the dynamics of a heavy particle that interacts with an ideal Bose gas of light 
particles we consider the Schrodinger equation 

id t %(x, yi ,...y N ) = H¥ t (x,y K ..y N ), (1) 

with 

A N N 

H:=-f+ pV(x) - £ A yk + J]w(x- y k ), (2) 
P k=\ k=\ 

for vectors in the Hilbert space of complex- valued, square-integrable (A^ + l)-particle wave 
functions on configuration space R 3(W+1) , i.e., 

K :=L 2 (R 3(/V+1) ,C). 

Here iel 3 represents the position of the tracer particle, and yi, . . .,y N 6 K 3 are the positions of 
the N gas particles. The Laplace operators w.r.t. x and y& are denoted by and A w , respectively. 
The potential W describes the interaction of the tracer particle with the gas, and the potential V 
describes an external force driving the tracer particle. The support of each gas particle tensor 
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component of ¥^=0 is given by a compact region A c R 3 the volume of which is denoted by |A|. 
The mean particle density is given by 

and is kept fixed. The ratio between the mass of the tracer particle and the one of the gas particle 
is chosen to be 2p. We use units such that h = 1 - note that V has different dimension than W. 

One might expect that if the tracer particle is well-localized initially, and because of its heavy 
mass, the variance of its position remains small. Therefore, its mean position should effectively 
obey Newton's second law, for a force generated by the external potential V and by the mean-field 
version of the potential W originating from an effective dynamics of the gas. It should be observed 
that, for large volume |A|, the influence of the tracer particle on the gas is not of leading order. Yet 
the back reaction of the perturbation of the gas on the tracer particle is relevant. In our study, we 
are interested in the (one-particle) wave function of the gas, as compared to a reference state </>f et) 
that obeys 

id t ^(y) = -A y 4 e %) (4) 

and describes the effective one-particle wave function of a system of gas particles in the absence 
of the tracer particle. The effective one-particle wave function of the system of gas particles inter- 
acting with the tracer particle is denoted by 

The function e, describes an "excited state" of the gas. 

For large |A| and p, the Schrodinger dynamics given by (1), (2) is expected to be well approxi- 
mated by a semi-classical dynamics described by the equations 



id t e t (y) = (-Ay + W(X t -y))e t (y) + W(X t -y)$ ef \y), (5) 
^ = -VV(X t )-VW*\^(X t )-2%VW*^^j(X t ). (6) 

Here ( h I, e R 3 denotes a classical trajectory of the tracer particle. Equation (5) is a Hartree-type 
equation, and e, is called Hartree (one-particle) wave function. In order to keep the excited state of 
the gas described by e t visible in the scaling considered in this paper, we choose the norm of <^ ref) 
such that the inhomogeneity in (5) is of order 0(1), i.e., 

lk£?IL = on). (7) 

Remark 1.1. In our case 0| ref) varies only little in a neighborhood of X t . Hence, it is possible to 
replace the inhomogeneous term in (5) by 

WiXt-y^Xx,). 

From now on, we refer to the time evolution generated by (1) as the microscopic dynamics and 
to the one generated by (5)-(6) as the macroscopic dynamics. The goal of this work is to quantify 
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the closeness of these two dynamics and to estimate the rate of convergence, as |A| ,p — » oo. 

Derivations of such mean-field equations from the microscopic Af-body Schrodinger evolution 
are usually carried out by making use of reduced density matrices and are based on hierarchies 
[Spo80, EY01]. In recent years, alternative methods have been developed to derive the Hartree 
equation from the microscopic dynamics. One approach, developed in [FKS09], relies on the 
Heisenberg picture and involves dispersive estimates and the counting of Feynman graphs. Another 
one was introduced in [RS09]. It exploits the dynamics of coherent states and is inspired by a 
semiclassical argument given by Hepp [Hep74], who initiated all these studies. In this paper we 
follow a different approach introduced in [Picl 1], which is based on counting the number of "bad" 
particles, i.e., particles that are not in the state given by the Hartree wave function. 

The Hartree wave function is of interest because it can be used to investigate physically interest- 
ing phenomena, such as quantum friction by emission of Cerenkov radiation [FGS1 1], more easily 
than by using the microscopic dynamics. Rigorous control of the microscopic time-evolution in 
terms of a macroscopic one is however a difficult problem. Our paper provides a first result in 
this direction for a system consisting of a tracer particle interacting with an ideal Bose gas. For 
the analysis of Cerenkov radiation, i.e., the deceleration of a particle with a speed higher than the 
speed of sound in the gas, interacting Bose gases, where the speed of sound is non-zero, are most 
interesting. The techniques to treat the ideal gas presented in this paper appear to be very robust 
and to allow for many generalizations. A mean field pair interaction of the gas particles, for ex- 
ample, can readily be introduced into our mothods, using estimates provided in [Picll]. Further 
generalizations to models with a thermodynamic scaling of interacting gases are presently under 
investigation. 

Acknowledgments: D.-A.D. gratefully acknowledges financial support from the post-doc pro- 
gram of the DAAD. J.F.'s visit at IAS is supported by 'The Fund For Math' and 'The Monell 
Foundation'. A.P. is supported by the NSF grant #DMS-0905988. Furthermore, D.-A.D. and P.P. 
would like to thank the Mathematisches Institut der LMU Miinchen and the Department of Math- 
ematics of UC Davis for their hospitality. 

2 Notation 

1 . The expectation value of an operator O w.r.t. the microscopic wave function % is denoted 
by 

<0> ( :=<T (> 0¥ ( >. 

2. |-| is the canonical norm on C d , for any dimension d; \\-\\ p is the norm on the Lebesgue space 
L p , < p < oo. For operators, O, on L 2 or *H, we denote by \\0\\ their operator norm. We 
also introduce the norm ||M|| p for matrix- valued functions, M, with matrix elements M,y. It 
is defined by 

\\ M \\p ■= Z INI ■ 

•J 

3. The momentum operator of the tracer particle is denoted by 

P ■= -N x . 
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Furthermore, for n e N, D (A") denotes the natural domain of the self-adjoint operator A" in 
L 2 (R 3 , C), and D (A 00 ) := n„ eN D (A"). Similarly, we denote by D(H ) and D{H n Q ), neN, the 
domains in *H of the self-adjoint operators 



A A 



H o-=-^-J]K and HI 



2 P ,= , 

respectively, and D(H£) := n neN D(tf£). 

4. The Fourier transform of a function 77 6 L 2 is denoted by 77. 

5. Given a vector 77 e L 2 with ||?7|| 2 = 1 we denote the orthogonal projection onto 77 by 

\n)(ri\. 

Furthermore, we use the notation 

q\ := l-pl (p v k y) (x,y u ...,y N ):= fj(y k ) J d 3 y k r]*(y k mx,yu ■ ■ .,y N ), I < k < N, 

(8) 

and 

q r } := 1 - ?r\ p ( r f) ■= H m # f) ) (iap i/2 ^ ref) i . (9) 



We note that, in q n k , the subscript k always stands for the /c-th tensor component, while, in 
q^, the subscript t always refers to time t. 

6. The convolution of two functions f,g on R 3 is denoted by (/ * g)(-) := j dy /(• - y)g(y). 

7. The statement "F e Bounds" refers to a continuous, monotone increasing function F : R + —> 

R + . 

8. If not specified otherwise the symbol C denotes a universal constant the value of which may 
change from one line to the next. 

3 Main Result 

Throughout this paper we assume that the density, p, of the gas fulfills 

p>l. 

In order to compare the solutions of the microscopic and macroscopic dynamics, we choose 
initial conditions that are compatible with each other: 

Definition 3.1. 
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(i) As initial conditions for the microscopic dynamics (1) we consider *P f |,=o = 'F®, where 



^\x, yi ,...,y N ) :=x m (x)Y\<P i0) (y k ), P { % = 1, 



(10) 



is given in terms of some unit vectors ^ (0) , </> (0) 6 L 2 (R 3 , C) with the properties 

X m eD(A°°), f'eDfA 00 ) 

such that: 

(a) The initial variance of the position x and the velocity 11 of the tracer particle fulfills 



,(0) 



+ 



'P-^PX™)) 
P 



^ 0) )<- x 



c_ 

-id 



for a fixed constant 6, with < 6 < 1 . 

(b) The support of (O) is contained in a compact region A c K 3 , and, as |A| 
initial wave function becomes flat, in the sense that 



0(0) 



= C|A| 



-1/2 



VftW <C\A\ ~ 5/6 . 



(11) 



oo, the 



(12) 



(ii) As initial conditions for the macroscopic dynamics (5)-(6) we choose 

£t\t=0 

Furthermore, we define 



= o, x t u = (x {0 \ x X {G) ) , = (x (0 \ ■ (13) 



to be the solution to (4) with 



(ref) . 



2/12)3 



(ref) 



^ ref) Uo = |A| i/z ^ u; . 



- I A 11/2^(0) 



(14) 



Remark 3.2. (i) Note that an example of a function ^ (0) satisfying (11) is a Gaussian wave packet 
with a variance in the position comprised between p~ y and p 7 ~\ for some < y < 1. By Heisen- 
berg's uncertainty principle, the variance in the momentum is then between p 1 y and p r , and, since 
the mass of the tracer particle is of order p, the variance of its velocity is between p _r and p 7 " 1 . 
Hence, 6 can be chosen to be min {y, 1 - y}. (ii) A product wave function like (10) is of course a 
very special initial condition. However, this condition can be relaxed, as pointed out in Remark 
3.6 below. 

In order to keep our analysis simple, we assume the potentials V, W to be smooth functions of 
compact support, i.e., 

V, WeC^(R\R). (15) 

It is a standard result that the equations of motion (1), (4), and (5)-(6), with initial conditions as 
chosen above, have smooth solutions. 
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Definition 3.3. We denote by 



¥ : R -» L 2 (R 3 ,C) ®L 2 (R 3 ,Cf N c "H, ? h> (16) 
the unique solution to the Schrodinger equation (1) with initial condition given by (10); by 

0(«f) : R->L 2 (R 3 ,C), t^^ ef \ 
the unique solution to Eq. (4) with initial data as in (14); and by 

£XXxX:R^L 2 (]R 3 ,C)xR 3 x]R 3 , t ^ (e t ,X t ,X t ) (17) 
the unique solution to equations (5)-(6) with initial data as in (13). 

Remark 3.4. Note that assumptions (15) and / ^ (0) , </> (0) e Z)(A°°), in Definition 3.1, are much 
stronger than necessary. As can be seen from the norms used in our proofs, the results presented 
below hold for a considarably more general class of potentials and initial wave functions. Finding 
optimal conditions is, however, not our aim in this paper. 

The main result in our paper is the following theorem: 
Theorem 3.5. Let 

P { r m) ■= q ( r f) v,( iM i%> <%i ) q ( : si \ p ( r m) ■= i*> <*i 

denote the density matrices of the gas excitations w.r.t. 0f ef) of the microscopic and the macro- 
scopic descriptions, respectively; for the definition ofq ( * ef \ see (9). There exist C\,C% e Bounds 
such that, for allt eR and sufficiently large p and | A|, the following estimates hold true: 

(i) 

|| p (micro) _ p (macro)|| < ^ ^1/2 +p -l/2 + |A| -l/3j (lg) 

<Ci(0(p"' min{1,fl + |Ar 1/3 ). (19) 

(ii) 

\\X t - (x) t || + \\X t - l^j || < C 2 (t){a 1 ' 2 +p- 1 ' 2 + |A|" 1/3 ) (20) 

< C 2 (t)(p-? mki[1 ' s] + |A|~ 1/3 ). (21) 

Here ao, defined in (25) below, reflects the dependence on the initial condition (10). 
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The operator tr xy ^_ yN stands for tracing out the degrees of freedom specified in the subscript. 
Estimate (i) in Theorem 3.5 quantifies how well the density matrix of the gas excitations w.r.t. 
the reference state 0| ret) is approximated by the effective density matrix \e t ) (e t \, while estimate (ii) 
quantifies how close the expected position and velocity of the tracer particle are to the classical 
ones. Hence, for large p and |A|, the microscopic dynamics and macroscopic dynamics yield arbi- 
trarily close predictions, and, for practical purposes, one may thus use the macroscopic equations 
to study the behavior of the system. 

Remark 3.6. (i) We remark that inequalities (18) and (20) hold for any initial conditions such that 

a = 0| pHco (l) + 0| AHoo (l), 

while inequalities (19) and (21) only hold for the initial conditions (10) satisfying (11) and (12). 
(ii) Theorem 3.5 can easily be generalized to systems of M > 1 interacting tracer particles, as 
discussed in Remark 4.6 below. 

4 Proof of Main Theorem 

The strategy of the proof is a two step procedure. First, we probe how well the gas particles in 
retain the product structure encoded in the initial wave function (10). In particular, we compare 
the microscopic dynamics to the effective dynamics generated by 

id t ip t (y) = (-A y + W((x) t -y))<p t (y) (22) 

for the initial value 

<P,\,=o = (O) - (23) 

Definition 4.1. We denote by 

<p :R—* L 2 (R 3 , C), t^(p t 
the unique solution to (22) with initial condition (23). 

As shown in [Picl 1] for pure Bose gases without a tracer particle, it is convenient to control the 
deviation of the gas wave function from a product wave function with the help of a Gronwall-type 
estimate of the form 

—a t <Ca t + —, where a, := ) { . (24) 

The quantity a t counts the relative number of tensor components in T^ not showing product struc- 
ture. This can be seen best by means of the identity (5) in Lemma 2.2 of [PiclO], i.e., 




where O denotes the symmetrized tensor product (see (8) for the definition of p^ and qf). Hence, 
a t corresponds to the expectation value of the ratio k/N between the number of particles, k, that 
are not in the state of the Hartree wave function and the total number of gas particles N. 
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The situation considered in this paper is more complicated because of the presence of the tracer 
particle, which couples to the gas and generates entanglement between its state and the state of the 
gas particles. As a consequence, the error estimates will not only depend on a, but also on the 
position and momentum distribution of the tracer particle wave function. To see this we consider 
the example of an initial wave function x <0) of the tracer particle formed by a superposition of 
two wave packets whose supports are separated by a distance of order one. In the worst case the 
mean position (x) t could then be somewhere in between the supports of these wave packets. In 
a situation like this, the effective interaction term W {{x) t - y) in (22) has nothing to do with the 
actual interaction given by Yuk=\ W(x - yt) in (2), and there is no reason to expect that a t stays 
small. 

Moreover, as discussed above, a t is the expected ratio k/N of "bad" gas particles w.r.t. to 
N. Yet, for the control of the dynamics of the tracer particle, this ratio will not be very relevant, 
because, due to the support of W, the tracer particle only sees 0{p) many gas particles at a time. 
In the worst case scenario however, even though k/N might be small, all the "bad" gas particles 
could actually be in the vicinity of the tracer particle. It is therefore important to know how many 
gas particles are "bad", as compared to a number of gas particles of 0(p). The latter amounts to 
estimating the quantity |A| (q^^j that gives the expected ratio k/p. 

We must therefore carefully adapt a, to our situation. It turns out that, among appropriate 
choices that make the desired estimates quite easy, the following one is convenient. 

Definition 4.2. We define 



,2,2 



«< ■■= j ((* - <^) 2 ); + ((^) ) t + (iM «) 2 + (im 2 («r# )/ (25) 

for t e R. 

Note that the function 1 h-» a t is smooth. A key estimate is the following lemma. 
Lemma 4.3. There are Q, , € Bounds such that, for all t > 0, the following estimate holds: 

ja t <CV(t)a t + C%\t)p-\ 

In a second step we then control the error made when replacing the mean position t i-» (x) t in 
(22), which fulfills the Ehrenfest equation 

^ (x) t = (-W(x) - X - V g W(x - y k )j , (26) 

by the classical trajectory t h-> X t obeying (6). Furthermore, in order to probe the excited modes 
of the gas, we need good control on how well the effective wave function of the gas, \A\ l/2 (p t , 
approximates 

For this second step, too, we will invoke a Gronwall-type estimate. We will consider the 
following expression. 
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Definition 4.4. We define 



t~ fit := fX t - {x) t f + ( d(Xt ~ t (X)t) ) + |(0f ef) + e) - Al 1/2 <p t 

for t 6 R. 

Note that the function 1 1-» is smooth. We will prove the following estimate. 



(27) 

2 



Lemma 4.5. There are C B , Cf 1 £ Bounds such that, for all t > 0, the following estimate holds: 



j f fi t < (f{t) (fir +fi r 2 ) + (fit) (^ t + a t + |A|- 1/3 ) . (28) 

This will complete the second step, and combination of Lemma 4.3 and Lemma 4.5 will com- 
plete the proof of our main theorem. 

Proof of Theorem 3.5. With the help of Gronwall's Lemma, Lemma 4.3 implies that there is some 
C a e Bounds such that 

a t <C a (t)(a +p- 1 ). (29) 
Moreover, the function fi t is smooth, and, by our choice of initial conditions 

fi,\,=o = 0. 

We assume further that 

a t \t=o = 0| p hoo(1) + 0|a Hoo (1), (30) 

i.e., that the right-hand side of (30) becomes arbitrarily small for sufficiently large parameters p 
and |A|. Hence, there exists a time T > such that, for all < t < T, the bound < fi t < 1 holds. 
The estimate in Lemma 4.5 then implies 

j t fi t < 2C^\t)fi t + Cf(t) ( Jc a (t) (oro + p" 1 ) + C a (t) (ao + p" 1 ) + |AP 1/3 ) , (31) 

so that Gronwall's Lemma guarantees the existence of some Cp 6 Bounds such that 

fi t <C p (t)(af+p- ll2 + \A\- 11 *). (32) 

Let T(p, |A|) be the supremum of all such times T, and let us assume that T(p, |A|) is uniformly 
bounded. Then, upon choosing p and |A| large enough in (32), we can arrange for fij(p\A\) < J- 
However, from Gronwall's Lemma it follows that (32) holds for some time t > Tip, |A|). This 
contradicts our assumption, and hence 

lim Tip, |A|) = oo. 

p,|A|— »oo 

We therefore conclude that (32) holds for arbitrarily large t > 0, provided p and |A| are large 
enough. This proves claim (ii) of Theorem Theorem 3.5. 
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Next, we prove claim (i). Similarly to [Picl 1], we start by inserting the identity 1 = pf + q 1 ^ 
on the left and right-hand side of \m t ) Q¥ t \, which yields 



(micro) (macro) 1 1 
Pt 'Pi 



< 



\qT f) tr x ^ yN \ x V t ) <^|^ ref) -|6,)<6 r ||| 
|A| qf' tr x , n _ yN [p\> \%) (%\ p» } q ( r f) - \e t ) (e t \ 
+2 |A| |^ ref) tr x , y2 ,..., yN [p*> \V t ) (%\ qf] q™ 

+iAi||?r^ >M [?ri^><^kr]^ ref) | 

In order to estimate (33) we need the preliminary bound 



(33) 
(34) 
(35) 



^ ef) |^)<^|^ ef) -|6,)<6,| 



(ref) 



wr- n - (<f>r + e t ) + (<r } + *)) (im i/2 <p t - (# f) + *) + ^ ref) + *)| ^ - 1*> ^ 



< 



| A |l/2 ^ _ + 6 J + 2 || |A| l/2 ^ _ ^(ref) + f J ^ + ||^(ref) ^ ^ ^(ref) _ fe) ^ 



</? 2 + 2A + ||pf ef) \e t ) /HI + 2 |^ ref) k t > <*| 



C e (0 2 . C e (t) 



+ 



|A| |A| 1/2 ' 



(36) 



where in the last two lines we are using (27) and Lemma A. 2. Next, we establish estimates on the 
terms (33),(34),(35): 

Term (33): By Fubini, one finds that 

(<Pt\ tr x ,y 2 ,..., yN [PF ( ) (%\] \<P,) = (%, pf%) = 1 - q*¥ t ) 

so that 

(33) = |||A| q ( r f) \<p t ) (<p t \ tr x ,y 2 ,..., yN [\%) (%\] \<p t ) (<pA q^ ~ k> 

(l - (%, ql'%)) [|A| q^ \to) (to\ q ( r f) " l*> <*l] + l*> <*l| 

(l - (%, tfv t )) [|A| q ( r i] \<p t ) (<p t \ <z! ref) - kr> + |(%, Ml 



< 



< 2 



'0l + 2Pt+ Wf + C ^ 



|A| |A| 1/2 J 
where we have used inequality (29) to get 

l-(%,qfW t ) 



(37) 



< 1 + — < 2 

|A| 



for |A| ,p » 1, and, furthermore, inequality (36), definition (25), and Lemma A. 2. 
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Term (34): 



(34) = 2|A| 



tr x , y2 ,...,y N [pf \%) (¥A qf] q\ 



ref) 



= 2|A| 



,(ref) 



(ref) + 6, 

|A| 1/2 



|A| 



1/2 







+ 


|A| 1/2 



S(ref) 



|A| 



1/2 



< 2|A| 

< 2V^(C e (0+A) 

where we have used definitions (25), (27) and Lemma A. 2. 

Term (35): 



+ e, 




(38) 



(35) 



|A| 



ef) 



< a„ 

where we have used definition (25). 

Collecting estimates (37), (38), (39), and using (29) as well as (32), we find that 

II (micro) (macro)|| fl 2 , to , C £ {t) Ceif) OC t 2 , ,-. / — , n \ , 

\\Pt ~Pt \\ < fi t +Wt + —TTT + 777777 + 77T C f W + 2 (C 6 (0 + /?*) + or, 



(39) 



< 



|A| 

c l( oK /2 +p- i/2 + 



|A| 1/2 |A| 
|Ar 1/3 ), 



for some Ci 6 Bounds, as well as 



Iix,-<.v),|| + ||x,-^ || < A 



< 



C /3 (t)(a l Q /2 +p- l/2 + 



|A|- 1/3 ). 



(40) 

(41) 
(42) 



where inequality (41) follows from the definition of the function [3, given in (27). The choice of 
initial conditions (10), satisfying conditions (11) and (12), ensures that 



C 

a t \ t=Q < — . 
P b 



(43) 



Hence, (43) together with inequality (40) proves claim (i), and, together with inequality (42), 
proves claim (ii). □ 

In the rest of this section we present proofs of Lemma 4.3 and Lemma 4.5 which provide the 
estimates on the time derivatives of the quantities a, and f3 t defined in (25) and (27). 

The aim in each of these estimates is to show that either the corresponding terms can be 
bounded in terms of a t and f3 t or that they are small if one of the parameters |A| or p is large. 
The main mechanisms exploited in our strategy are: 
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• Expansion of differences for smooth functions / of the form 

f(x) -f((x) t ) = R f (x, (x) t ) ■ (x - (x) t ) 
where R f denotes Taylor's remainder term. This is meant to take advantage of the estimate 

\\(x-(x),)¥ t \\ 2 < V^- 



• Rearrangement of the arguments of scalar products, if necessary by inserting the identity 
1 = pj' + q^, such that the operator acts directly on TV This is meant to take advantage 
of the estimate 



IM 



1/2 - 



• Whenever the Hartree wave function ip t is integrated against a function / with support con- 
tained in a volume of 0(1) one gains a factor |A| _1/2 , by Lemma A. 1 . Note that, for example, 



nfPt 



P\'f(x-yi)PT = (f*\<Ptn(x) P \, 

and therefore, thanks to the argument above, one gains a factor |A| _1 . 
Proof of Lemma 4.3. Since a t is a smooth function of t, we can estimate its derivative by 



dt 



d 



+ 



d_i(p~ (P) t 
dt 



+ |A| 



+ |A| 2 



jM)\ 



(44) 

(45) 

(46) 
(47) 



Denoting commutators by [•, •] and anti-commutators by {•, •}, we shall use the following auxiliary 
computation in the estimates of the individual terms (45)-(47): Let A be an arbitary self-adjoint 
operator; then the estimate 



|((A-<A>^)J = \i([HAA-(A) t f]) + ^ t (A-(A) t )^ 



d 



i ({[H,A - (A) t ] ,A - (A) t }) t -2 (A- {A) t ) t - (A\ 



< 



\({[H,A],A-(A) t }) t \ 
2\([H,A](A-(A) t )) t \ 



(48) 



holds true, supposing the expressions on the right-hand side are well defined (recall that is nor- 
malized). 
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term (44): Using definition (25) we estimate 



(44) = 



< 2 



P ~ <p) t 



(x - (x) t ) 



< 2 

< Ca, 

term (45): With the help of (48) we get 



P ~ (P)t 



{{x-{x),f) t 



(45) = 

< 2 
= 2 

< 2 



d_ l(p-(p) t 
clt 



H, — 



P ~ (P)t 
P 



pV(x) + NW(x - yi), - 
P 



P ~ (P) t 
P 



VV(jc) 



P-(p) t \ 



+2\A\ 



yW(x- yi ) 



P~{p), s 



We expand VV(jc) according to 

VV(jc) = W((x) t ) + R vv (x, (x) t ) (x - (x\), 
where R vv denotes Taylor's remainder term, and, using (25), we obtain 

l P~{p)t 



(50) < 2 



W«*>,) 



+ 2 



(R vv (x,(x) t )(x-(x) t )) 



P~{p), S 



< 2\\R vv \\j{x-(x) t )V t \\ 2 



P ~ (P)t 



(49) 



(50) 
(51) 



(52) 



The estimate on term (51), which depends on W, is more involved. It is convenient to split it as 
follows, inserting the identity 1 = p^' + q^', 



(51) < 2|A| 
+2|A| 
+2|A| 



P 



(pl'+q^VWix-yJq* 1 



P-(p) t ] 



qV'VWix-y^-?—^ 



(53) 
(54) 
(55) 
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In order to treat the term (53), we expand VW according 

VW(x -y) = VW((x\ -y) + R vw (x, y, (x) t ) (x - (x) t ), 
where R vw denotes Taylor's remainder term, and find 



(53) = 2|A| 
= 2|A| 



VW * \ip t \\x)p\ 



I 



2,....,* ES 

dyVW((x) t -y)\<p t \ 2 <y)-(p* 



P 



P-(P), ] 
P I 



+2\A\ 



J dy [R vw (x,y, (x) t ) (x - (x\)] \<p t \ 2 (y)pf 



P-(p) t ] 
P l 



To estimate term (56) we use definition (25) to get an auxiliary bound 



P~{p)t S 
P I 



P ~ (P) t ^ l^p, p-(p) t y 



ft 



P ~ kp)t 
P i 



p - (p) t 



< 



< 



|A| 
a, 



p - (p)t 
p 



|A| 



1/2 ! 



which implies 



1/ 



(56) < 2|A| dyVW((x) t -y)\<p t \ z (y) 



|A| 

2 



1/2 



< 4 iai (HvwiL + www,) [ ikvi: + ii^ - e i&t rv — 



|A| 



1/2 



— C Cprop(0 



2 a t 



|A| 



1/2 



Here we have inserted the identity 

(p t = e (p + ip t - e qr 
and used Lemma A.l, which implies that 

IkVT + Ik - < <wt) . 

To control term (57) we make use of definition (25) and Lemma A.l again and find that 



(57) < 2|A| 



< 2|A| 



dyR^(x,y,(x) t )\<p t \\y) 



sup I 

VxeR 3 J 



(x-(x) t )%\\ 2 



n . P~(P)t xv 
P\ ~t 



dyR vw (.x,y,(x) t )\<p t \ 2 (y) 



((x-(x) t ) 2 \ 



P~{p)t 



P 



< 4|A| sup (||^V,sz)|L + ||^ v "Usz)||JI|r^||; + ||^-e-y 



WW 



— C Cp rG p(?) 
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Estimates (58) and (59) imply the upper bound 



(53) = (56) + (57) 

The term (55) is estimated similarly: 

c, P - (P)t ] 



(60) 



(55) = 2|A| 



q^W{x-y x )p\ 



< 2\K\\\q^ t \\pW( < x-y x )p\ 



P ~ (P) t 



P 



< 2 



|A| |AP 1/2 ^j\A\ Jfi) t pW(x- yi )p»\ 



P ~ (P) t 



< C|A||A|- 1/2 yf^ t sup(\\VW(x--)\L + \\VW(x--)\\ 2 )(\\e iAt <p {0) \l+ 

*sR 3 

< C |A| |AP 1/2 V^C prop (f) |A|" 1/2 

< CC pmv (t)a t , (61) 
where we have used definition (25), Lemma A.l and the following straightforward estimate 

pW(x- yi )pff = sup tap? [VWix-yOfp^'x) 

IMI 2 =i 

= sup l X ,pl' f dy[VW{x-y)} 2 \ip t \ 2 (y)p^x) 
< sup f dy [VW(x-y)] 2 \if t \ 2 (y). 

jteR 3 J 

Our estimate on term (54) is more subtle. Due to symmetry in the gas degrees of freedom, defini- 
tion (25) of ar f , and (3) one finds that 



(54) = 2|A| 

= 2|A| 

< 2p- 1 

< 2p~ l 



VW(x- yi )q\ 



n P ~ (P)t 



N 

J]ql'VW(x-y k )% 



k=l 

N 



P ~ (P)t 



1/2 



£ VW(x - y k )ql'q^ ■ VW(x - yj) 
kj=\ 

< P~ 2 {Yj V W (* ~ ytekJj • VW(x - yj )j + a t . 



(62) 
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In the next step, the sum is split into a sum over diagonal terms and one over cross terms, i.e., 

N 



p~ 2 1 £ VW(x - y k )q*tf -VW(x- y 7 )| 

= p- 2 N(vW(x-y 1 )q' p 1 '-VW(x-y l )) t (63) 
+p~ 2 N(N - 1) (VW(x - yi)q?q$ ■ W(jc - y 2 )) . (64) 

The diagonal term (63) can be bounded by 

(63) < p- 2 N\\VW(x- yi m^\\\\VW(x-yiyift\\2 
= p- 2 N((VW(x- yi )) 2 ) t 

= p- 2 N(p*'(VW(x- yi )) 2 p*') t (65) 
+2p- 2 N 3 (pf (VW(x - yO) 2 <) ? (66) 
ViV(^'(Viy(x-yi)) 2 ? f) ( , (67) 

where we have again inserted the identity 1 = p^ + q*. Using definition (25) of a t , as well as 
Lemma A. 1 , we compute 

(65) < Cp- 2 N sup (\\VW(x - + - 0lQ(||e ,7 V O) ||L + |k - e ! ' A y 0) " 2 



-prop 

and 



< CC prop (0V 1 , (68) 



(66) = 2p' 2 N3(p v l '(VW(x-y 1 )) 2 q' p l ') t 

< 2p~ 2 N \\(VW(x - y,)) 2 pf%\\ 2 \\qf%\\ 2 

< Cp~ 2 N sup (||(VW(x - -)) 2 L + ¥ VW ( X ~ ')) 2 U x 

jceR 3 

x(||^°1L + lk-^ (0) || 2 )iAr 1/2 v^ 

< c P - 2 ivc prop a)iAr 1/2 iAr 1/2 

< cp _1 c prop (oV^ 

< CCp.opWfo'r+p" 2 ), (69) 



and finally 



(67) = p- 2 N(q^(yW(x- yi )) 2 q v ') t 

< p- 2 N\A\~ 1 ' 2 V^||(VW) 2 |L|Ar 1/2 ^ t 

< Cp~ l a t . (70) 
To estimate the cross terms in (62) we again insert the identity 1 = p^' + q[' on the right and the 
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identity 1 = p 2 + q 2 on the left of the expectation value (64). This yields 

(64) = p~ 2 N(N - 1) + q\')q^VW{x - y{)VW(x - y 2 )q*(p% + q*)\ 
< p~ 2 N(N - 1) \(p"q%VW(x - yi)VW(x - y 2 )q\'p 2 H ) t 
+2p- 2 N(N - 1) \(p?q%VW(x - yi)VW(x - y 2 )q\ t q%) t 
+p- 2 N(N - 1) \(q?q*VW{x - yi)VW(x - y 2 )q\'q%) t ■ 
Each term on the right-hand side can be estimated using definition (25) and Lemma A.l: 

(71) < p~ 2 N 2 \(q v 2 'p v l 'VW(x-y,)VW(x-y 2 )p v 2 'q v 1 ') t 

< p- 2 N 2 \\qf%\\ 2 \\VW(x-yi)pf\\ \\VW(x-y 2 )p%\\ \\qt%\\ 2 



< Cp~ 2 N 2 |A|" 1/2 x 



x 



sup (|| W(* - OIL + \\*W(x - -)|| 2 ) i\\e m <f ] \ 



+ \\<pt - 



x 



x|Af 



-1/2 



a t , 



and 



(72) < 2p- 2 N 2 \(q 2 'p v <VW(x-y l )VW(x-y 2 )q v l 'q 2 °') t 

< 2p~ 2 N 2 \\q^%\\ 2 \\pfVW(x - yi )\\ HVWIL 11^9^4 



< CpN 2 \A[ 



-1/2 



a, x 



x sup (\\vw(x - OIL + \WW(x - oil 2 ) ( |k' A y 0) L + h 



finally 



xiap 1 



(73) < p- 2 N 2 \(q\ t q^W(x-y l )VW(x-y 2 )q\ t q%) t 

< p~ 2 N 2 ||<^^|| 2 ||VW||i \\q^q^t\\ 2 

< Cp- 2 N 2 \A\~ l ^ t \A\~ l 

< Ca t . 



(71) 
(72) 
(73) 



(74) 



(75) 



(76) 



Combination of estimates (52), (60)-(62), (68)-(70), (74)-(76) implies the existence of a C Var (p) e 
Bounds with 

(45)<C Yaiip) (t)(a t +p- 1 ). (77) 

With the help of definition (25) of the function a, and Lemma A. 1 we treat the two remaining terms 
as follows. 
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term (46): Recall that, by (22), the equation of motion 

j/ 1 ' = 'l/ 1 ' = ~^~ Ayi + w ^-y^Pi] = i[-A» + WHxh-yiM?] 

holds. Using the expansion 



W(x - y) = W((x) t -y) + R w (x,y, <*>,) (x - (x) t ), 
where R w (x,y) denotes Taylor's remainder, we find 

(46) = |A| 

= |A| \([W(x - yi ) - Wdx), - y{), <])J 

= 2 |A| |3 ((pf + 9 f ) (W(x - yi ) - W((x) t - yO) q") t 

= 2 |A| |3 (pf (W(x - yi) - W«A - yO) 9 ?) f 

= 2|A||g(pf/? w (x,y 1 ,<x> f ) (x-(x) t )q (p l ') t 

< 2 |A| \\(x - (x) t ) %\\ 2 \\R w (x,y u <x>,)pf || ||<%|| 2 

< C | A| yfc sup (||i? w (x, z)\l + \\R w (x, ■, z)\\ 2 ) (\\e iAt 

«€]R3 

< C C pm p(t)a t , 

where, in the step from (79) to (80), we have used that 

q\'(W{x-y{)-W{{x) t -y{))q\' 
is self-adjoint, so that its expectation value has a vanishing imaginary part 



iAti(0) 



|A|- 



term (47): 



(47) = |Af 



= |A| 2 \([W(x - yO - W((x) t - y0 + W(x - y 2 ) - W((x) t - y 2 ), q*q* 

= 4 |A| 2 |g (ql' (W(x - yi ) - W((x) t - y0) q*q*)\ 

= 4 |A| 2 |3 ((pf + ^ (W(x - y0 - W«x>, - yd) q^q p <) t 

= 4 |A| 2 |g (pf ^ (W(x - yi) - W((x) t - y0) 

< 4|A| 2 (|(pf^'W(x-y 1 )«)| + |(p^2 w «^|-yi)«W>|) 

< 4 |A| 2 ||^4 (||W(x - yi )pt'\\ + \\W((x) t - yOpfH) hTql'%1 

< C |A| |AP 1/2 sup (\\W(x - OIL + W(x - -)|| 2 ) x 

+ ||^-e'- A y o) || 2 )iAr 



x(||e' A y 0) ||, 

< cc prap (0ff,. 
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Note that (82) holds, because every other operator inside the expectation value commutes with qt' . 
In the step from (83) to (84) we have used that the operator 



.ft >pt 



q v 'q v 2 ' (W(x - yi ) - W((x) t - yd) q*'q 
is self-adjoint, and hence the imaginary part of its expectation value vanishes. 

Upon collecting estimates (49), (77), (81), (85), we find C^, e Bounds such that 

d 



-a t <C^\t)a t + Cf{t)p-\ 
at 



which concludes the proof. 



□ 



Finally we prove the bound on the time derivative of fi t . The aim is to show that either the 
corresponding terms can be estimated in terms of the quantities a, and/or /3 t , or else that these 
terms are small if one of the parameters |A| or p is large. 

In addition to the arguments we have used in our proof of Lemma 4.3, we exploit the following 
facts: 



• The crucial cancellation 



pfVWix - y 1 )p v 1 ' - VW * \<p t \\x)p^ = 0. 



which will take place in (94) below. This cancellation determines the structure of the macro- 
scopic equations of motion. 

• The uniform spreading of the wave function $ ef , which is controlled by the propagation 
estimate (ii) in Lemma A. 1 . With its help one can show that the convolution 

VW * \(p ( r f) \ 2 (x) 

is small for large |A|. 

• The propagation estimates on the wave function e, of the gas excitations w.r.t. 0| ref) , which 
are provided by Lemma A. 2. 



Proof of Lemma 4.5. Recall that t h-> j3 t is smooth, so that 



-r (**-<*>,) 
at 



+ 



+ 



dt 2 



(X t - (x) t ) 



1/2 



(86) 
(87) 

(88) 



We shall estimate the terms on the right-hand side individually. 



Term (86): By definition (27) of the function /?, we immediately get 



(86) = 



— (X, - (x) t ) 
dt 
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Term (87): Using the equations of motion (6) and (26), we find 
d 2 



(87) = 

< 2 



dt 2 



(X t -(x) t ) 



|w(X t ) - VV(jc) + VW * \e t \ 2 {X t ) + 221 W * (<^ ref) e,) (X t ) _ 7^ V JE W ( x ~ 



Using the expansion 

VV(x) = VV(X t ) + R vv (x, X t ) (x - X t ), 
where R vv denotes Taylor's remainder, we obtain 

\(VV(X t )-VV(x)) t \ < \(R vv (x,X t )-(x-(x) t + (x\-X t )) t 



< \\R 



(x-(x) t )%\\ 2 + \\R™\\Kx) t -X t 



< Cyfa t + C(3 t , 
where we have used (25) and (27). Furthermore, 

V W * \e t \ 2 (X t ) + 2%VW * (^ ref) e f ) (X t ) - i V Jg W(x - y k )j 

= \(VW * |0| ref) + 6,| 2 (X f ) - VW * \(P ( r f) \\X t ) - |A| VW(x - yi )) f 
< |(VW * |# f) + 6,| 2 (X f ) - |A| VW * \<p t \\X t )) t 
+ |A| |(W * \cp,\ 2 (X t ) - VW * \<p t \\x)\ 
(vW*\cp t \ 2 (x)-VW(x- yi )) t 



(89) 



+ 



IM 

|vw*i^! ref) ta)|. 

All these terms are estimated below. In the step from (90) to (91), we have used the identity 

U(ref) , _|2 i^(ref)|2 



(90) 

(91) 
(92) 
(93) 
(94) 
(95) 



+ e,r - 



\e t \ 2 + 23ty<%, 



and symmetry in the gas degrees of freedom to replace 



1 N 

-VYw(x-y k ) 



by 

\A\VW(x- yi ). 
To estimate the term (92), we use the identity 



+ e t \ 2 - |A| |^| 2 = % 



(^ + e t + \A\ 1 ' 2 < Pt )(^ n + e ! -\A\V 2 v ! ) 
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which, together with definition (27) and Lemma A.l, implies that 
dzVW(X t - z)2t 



(92) = 



< 



(<f> { r%) + e t {z) + |A| 1/2 <p t (zj) (cP ( r f \z) + e t (z) - |A| 1/2 <p t (zj) 
VW(X t - ■)% (^ ref) (-) + e t (-) - |A| 1/2 (p t (.) + 2 |A| 1/2 ^(-)) 



(96) 



x||0r+^-iM i/2 ^i 

< C(HW|L + ||W|| 2 )x 



x (||^ ef) + e t - iAr 2 ^n 2 + 2 Ar 2 || e 'V 0) L + 2 > A i 1/2 Ik - ^ (0) ll 2 ) x < 97 ) 



xll^ + ^-IAI 1 / 2 ^ 
< c(/? 2 + C prop (f)A), 

where, in the step from (96) to (97), we have used the identity 

-TO) , „ JAl j(0) 

Next, by expanding according to 

W(x - y) = VW(X, - y) + fl vw W,X,) (X, - x), 

one gets 
(93) = |A| 
= IM 
< IM 



(98) 



j dy [VW(X t - y) - V W(x - y)] N 2 (y) 
I^J dyR vw {x,y,X t )\<p t \ 2 (y) (X t - (x\ + (x) t - x) 



dyR vw (x,y,X t )\(pr(y)(X t -(x) t )] 



< |A| sup 



+ IM 



dyR vw (x,y,X t )\<p t \\y)((x) t -x)) 



dyR vw (x,y,X t )\tp t m 



[\X t -(x) t \ + \\((x) t -x)V t \\ 2 ] 



< C|A| sup (\\R vw (x,;z)\\ + \\R^(x,;z)\\ I )(|k w m: + Ik 



_/Ar .(0)1 



,*A? ,(()) 



< C|A| C prop (0 2 |Ar 1 ( y 8,+ v^) 

< cc prop (o 2 (A+ V^), 

where we have used definitions (25) and (27), as well as Lemma A.l. 
Moreover, by inserting the identity 1 = + q^' and noting that 

p?VW(x-yi)p'* = VW* \<p t \ 2 (x)p*', 



(99) 
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we find that 



(94) = |A||((pf + q'*)[VW*\<p t ?(x)-VW(x-y 1 )](p? + q[')) t 
< 2 |A| |(< [VW * \<p t \ 2 (x) - VW(x - yi )] pf)J 



+ |A| 



< 



[VW * \<p t \ 2 (x) - VW(x - yi j\ q*) f 

2 |A| ||<%|| 2 (||VW * |^| 2 (jc)|| + pW(x - yO^'H) 
+2|A|||^Y ( || 2 ||VW|L|| 9 fT f || 2 



< ciai iap i/2 v^oivwil + iivwiii)(||e' A y u) L + Ik' - ^ (0) 

+2]A] |AP 1/2 V^I|VW|L|AP 1/2 

< CC prop (0 yfa t + Ca t . 

Finally, by integrating by parts, one arrives at the estimate 

(95) = |W*|<^ ref) l 2 (^)| 

< 211^11^11 NHL, 

which, with the help of Lemma A.l, yields 



(100) 



(95)<2C 2 f ||W%|A| 



-1/3 



(101) 



The estimates (89), (98), (99), (100), and (101) imply the existence of a function C ve i e Bounds 
such that 

(87) < C vel (0 (p t +tf + + a, + |AP 1/3 ) . (102) 
Term (88): With the help of the equations of motion (5) and (22) we find that 



(88) 



d 



1/2 



< 



23 ((-A v + W(X t - •)) (^ ref) + e) ~ ("A, + W((x) t - •)) A l/2 <p t , («^ ef) + e) - A i/2 <p t ) 



1/2 



1/2 



= 23 ((-A y + W(X t - •)) (& (ref) + * - A 17 2 <Pt) , <PT f) + e t - A 1 
+23 ((W(X t -■)- W((x\ - •)) A 1 V (0| ref) + e t ) - A 1/2 ^) 
Because the operator 

-A y + W(X t -y) 
is self-adjoint, the term in (103) is zero, so that 



(103) 



1 1/2 



(88) < |23((W(Z t --)-W«x> f --))A 1/ V# f) + ^-A 1/ V}| 

< 2 ||(W(X, - •) - W((x) t - •)) A^ 2 ^|| 2 /2 \\4f* + e t - A 1 *^ 

< 2 |A| 1/2 \\(W(X t -■)- W((x\ - •)) cpt\\ 2 + 2A- 



(104) 
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We expand W according to 

W((x) t -y) = W(X t -y) + R w ((x) t ,y,X t ) ((x) t - X t ), 
and use estimate (88) to obtain 

(104) = 2\A\ 1/2 \\R w ((x) t ,;X t )cp t (X t -(x) t )l + 2/3 t 



< C|A|^A sup (ii^c^, -,^>|L + ll^ w c^, -^)|| 2 ) (||^ A v (0) IL + 11^ - ^ A ^ (0) ll 2 ) + ^ 

< C(l + c prop (0)A, (105) 

where, once again, we have used (27) and Lemma A.l. 
Given estimates (102) and (105), one infers that 

j { Pt < Cf(t) (A +0}) + Cf\t) ( yfc t + a t + |AP 1/3 ) , 

for Cp\ C^ 2) 6 Bounds. This concludes the proof of Lemma 4.5. □ 

Remark 4.6. Note that our analysis easily generalizes to systems of M > 1 tracer particles with a 
microscopic Hamiltonian of, for example, the following form 



M 



(106) 



where we denote by xt, x%, . . . ,x M and p\, p 2 , ■ ■ . ,Pm the positions and momenta of the M tracer 
particles, and by / a regular pair potential. To prove a result analogous to Theorem 3.5, comparing 
the microscopic dynamics generated by (106) to the macroscopic dynamics, 



d ( ^ 
l dt £t(y) = 



M 



M 



■A y + J] w ( x kt ~ y) S(y) + J] W(X Kt - y)$°%), (107) 
k=l ) k=i 

df _ -W(X Kt ) - ^ VI(X K1 - Xj, t ) - VW * \e t \ 2 {X Kt ) - 2%VW * («^ ref) 6 f ) (X,, f ),(108) 



U Xfy 



for all 1 < k < M, the following natural adaptations are needed: 

1. As in Definition 3.1, we need to assume the initial wave functions 

*F (xi , X2, . . . , xm, yi,y%,---, yN) 

to be given by a product of a localized wave packet for the tracer particles - compare to x {0) 
in (1 1) - and a product wave function nf=i <P {G) iyk) for the gas particles. 

2. The appropriate intermediate dynamics of the effective wave function of a gas particle (22) 
is of the form 



id,<p t (y) 



M 



-A, + ^f((4-j) 



k=\ 



<Pt(y)- 
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3. The quantity a t defined in (25) must be replaced by 



M 



' k=\ k=\ 



Pk ~ (Pk) t 



+ 



(iai(^)J 2 +(iai 2 («)J 2 



4. The quantity J3 t defined in (27) must be replaced by 



|A|"V 



Taking advantage of the commutation relations 

[x fe , py] =0, V£ * j, 

Lemma 4.3 and Lemma 4.5 can be proven along the same lines as demonstrated for a single tracer 
particle. The only new terms are ones depending on the pair potential /, and they are of the form 



VI(x k - xj), 



Pk ~ (Pk) t 
P 



(109) 



Using an expansion of V/ of the form 

VI(x k - xj) = V/ (<**>, - (xj)j + R v k \x k , x h (x k ) t ) (x k - {x k ) t ) + R]\{x k ) t , x Jt (x } - , 

where Rj' denote Taylor's remainder terms, one can treat (109) in the same way as the term (50) 
in the case of only one tracer particle. 

A Propagation Estimates 

Lemma A.l. There are C prop , C re f € Bounds such that: 

(i) The solution 1 i-» <p t to (22) with initial value (23) fulfills 

py o) L + Ik - ^ (0) || 2 < c prop (o im- i/2 

for all times t > 0. 

(ii) The solution 1 1-» <^ ret) to (4) with initial value (14) fulfills 

II^IL * c «f' |K ef) L < C ref IM" 1/3 • (HO) 



Proof. 
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(i) Because of (23) and (12) one immediately gets the estimate 



\e iAt <f> {0) \ 



< 



g /A^(0) = 0(0) 



C|A| 



-1/2 



(111) 



for all t e R. Moreover, any solution to (22) fulfills the integral equation 

<p t = e iA y 0) - / f ds e iA ^ [W ((x) s - •)] <p s . 
Jo 

Hence, we infer the estimates 

f dse iA «- s) [W((x) s --)]<p s 
Jo 

< fds [W({x) s --)]U-e iAs ^)\ 
Jo 1 



Up t -e™<t>' 



i&A0)\ 



+ fds\\[W((x) s --)]e iAs <t> (0) \ 
Jo 



< \\w\i 



fds\\<p,-^><>™ 
Jo 



e l&s rT + t\\W\\ 2 sup e l&s cf> 



se[0,t] 



JAsAO) 



Inequality (111) and Gronwall's Lemma ensure the existence of some C prop e Bounds such 
that 

||^-e'V 0) || 2 <C prop (0|Ar 1/2 . (112) 
Estimates (111) and (112) prove our claim. 

(ii) Equation (4) and (14), together with the estimate (111), implies that 



Ik! 



(ref)l 



| A |i/2|L;Ayo)|| < c 



Similarly, with (12), one finds that 

||V^ ref) || = |A| 1/2 || e ' A 'V0 (O) || <|A| 1/2 £ ^'V^» 



IM 



1/2 



<C\A\ 



-1/3 



□ 



Lemma A.2. Let 1 i-» e, be the solution to (5) with initial data e t \t=o = 0. There are C e , C e e Bounds 
such that, for all times t > 0, the following estimates hold: 

(i) \\e t \\ 2 < C e (0. 



/ • • i (ret I - 

(n) M e,\\ < 



2 - |A| 



1/2- 



Proof. 
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(i) Since the homogeneous part of (5) is self-adjoint we may infer control over the norm of a 
solution t i-» e, from the inhomogeneity according to 



Jjj 

Jo 



W(X, - W^Lds 



< t \\W\\ 2 sup H^IL 

se[0,t] 

< tC C re f 
=: CM- 

(ii) Using the equations of motion (4) and (5), as well as initial condition (13), a direct compu- 
tation yields 



-i I w(x, - ■)£,) - i ( w(z s - o^r 



1/2 ' 



IM 



ds 



< 



< 



|A| ie[0,f] 



2 sup | k5 || 2 + -^sup |K'- e/) ||LllW% 

se[Q,t] \A\ se[0,t] 



-C re f C C e (f) + , TjCj-gf c 



|A| l/2 -r ,v , |A|1/2 

C 6 (Q 
' IM 1/2 ' 

which holds because of (1 10) in Lemma A.l. 



□ 
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